Localization of water monomers inside ice-like clusters. 
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Abstract 

On the basis of the experimental data we suggest that water monomers could be trapped 
in channels running through ice-like clusters in water. Our argument relies on a simple model 
that describes the motion of a dipole particle inside a channel in the presence of an electric 
field with linear gradient. The model admits of both finite and infinite regimes of motion so 
that the finite one could correspond to the particle being confined to a channel. 



1 Introduction: Water monomers trapped in channels of 
ice- like structures. 

It was W.C. Rontgen, [T] who suggested the presence of ice- like structural fragments in liquid 
water, and put forward the hypothesis of two fluid structure of water in equilibrium. Later, 
O.Ya. Samoilov, [5], and L.Pauling, [3], see also [1], extended the Rontgen model by introducing 
the concept of clathrate, low density and high density water. The idea has obtained experimental 
support that provides evidence for the existence of two kinds of water. Thus, using the technique 
of four-photon spectroscopy we have observed, for the first time to the best of our knowledge, the 
rotational resonances in water due to monomers H20, H2O2 and OH [5], [Q. It is important 
that the mobility of monomers could explain the high permeation of water channels, in biological 
membranes, [7]. As to the ice- like structural fragments, recently, Nielsson et al,|8], using x-ray 
spectroscopy, have found experimental evidence in favor of their existence. They visualized them 
as cluster structures, of several tens of A and hexagonal ice symmetry Ih. 

Considering water as a mixture of ice-like clusters and monomers, makes for better understand- 
ing the physics of water. Nevertheless, there are still questions about the nature of coexistence of 
the monomers and the ice-like clusters. Specifically, the penetration of the H2O monomers through 
channels of the ice-like clusters needs exploring. The key to the problem is provided by the study 
of beams of fast particles in crystalline media, effected many years ago for the needs of nuclear 
physics, [9j. The main point is that the velocity of a particle is greatly enhanced if the beam is 
directed along a crystalline axis providing a kind of channel. Thus, one can claim the channeling 
effect, which appears to be of a significance reaching outside nuclear physics. In fact the problem is 
the old one. Years ago, Ya.L Frenkel, [3], considered the motion of molecules and ions in channels 
inside crystalline structures, see also pH] . 



The nh.'se.rva.tinn nf rotatinnal mnrles nf water mnnnmers siiP'trests the existence nf cavities that 



bonds in water, that is 1 — 2 ps according to paper [TT]. The channels inside a cluster with 
hexagonal ice Ih structure, of diameter 5.7 A, could serve such cavities. It is worthwhile to notice 
that similar equipotential cavities can be provided through the freezing of water in a cryogenic 
matrix formed by argon, [T^], and carbon nanotubes of diameter 14 A, [T3]. It is reported. [T5] . 
that molecules of water preserve their mobility in carbon nanotubes down to temperatures ~ 8K. 

The main point about the dynamics of water monomers is that they could remain trapped in 
crystalline channels instead of forming complexes of hydrogen bonds and quitting the channels. In 
this paper we are considering a simple qualitative model that could accommodate the phenomenon. 
Our arguments essentially rely on the interplay of translational motion of a monomer and its 
rotational dynamics caused by the dipole moment of water molecule . 




Figure 1: Potential U = p E{x) cos{(p) 



2 A quasi-classical model for the motion of monomer 

We shall draw a qualitative picture of the motion of a monomer as the dipole in the following 
potential U — pE{x) cos{ip). Potential surface is presented in FigH]. To that end we shall consider 
the latter as a two-dimensional rotator moving in direction of axis-cc in an external electric field, 
E(x), which mimics constraints imposed by the channel. Thus, its configuration is specified by 
angle ip describing its rotation and coordinate x its position in axis- a;. Its dynamics is described by 
angular momentum L and momentum P. It is important to choose the right characteristic scales 
of the system. We shall take: 

(1) time scale t = 10^^"* sec corresponding to rotations of water molecules; 

(2) spatial scale r = 3 • 10~^ cm, close to the size of a water molecule; 

(3) mass scale 3.2 • 10^^'^ corresponding to a water molecule; 

(4) dipole moment ID = IQ-^^ CGS; 

(5) electric field 10 kV/cm. 
Hence we infer: 

(6) velocity and momentum scales 3 • 10^ cm/ sec, lO"'^^ gr cm sec~^ , respectfully; 

(7) moment of inertia « 3 • 10~^^grcm^. 

We shall assume that the external electric field, E{x), is linear in x 



E{x) =Ax + W 



The scales are conducive to the use of quasi-classical approximation and numerical simulation. In 
fact, for the above characteristic scales we have the de Broglie wave length 



10"" cm. 



Assuming the size of a channel C ~ 10 A, we get the parameter of quasi-classical expansion 



A 
C 



^ = - « 10" 



that is sufhciently small. The usual constraint imposed on the interaction potential, |14j . 



dU(x) 



dx 



> 



Aa;2 



which means that the potential is smooth enough at the de Broglie length, is valid for the interaction 
between the dipole moment and the external field 

U{x,ip)—p- E{x) cos{(p) 

for E{x) is linear in x. The similar requirement for the rotational motion is also verified. We have 
the dimensional expansion parameter 



= y w 3 • 10 



L = /(^ « 3 • 10 erg ■ sec 



It should be noted that the requirements indicated above are not satisfied at return points 
where the size of de Broglie length rapidly changes. For example, this is the case even of the 
harmonic oscillator. Thus, the conventional quasi-classical approximation breaks down, [14j . At 
these points one can change the x-representation for the p-representation, which, from the formal 
mathematical point of view, amounts to the Fourier transform. But the advanced theory, |15) , |16) , 
of the quasi-classical approximation claims that one can still use the classical equations of motion 
in conjunction with the Bohr-Sommerfeld quantization condition 



pdq = 2 7rnh, where n is an integer, 
provided the gradient of the Hamiltonian is not degenerate 

^0, (1) 



dx J \ dp 

where x, p are coordinates and momentum of the system. It should be noted that we are not 
constructing the wave function at turning points, but only find the quasi-classical trajectory by 
means of the classical Hamiltonian equations. It is the situation of the 'old' quantum theory by 
Niels Bohr. 

The arguments given above enable us to describe the dynamics of rotator within the framework 
of Lagrangian mechanics. The Lagrangian function reads 

C ^ ^-ip^ + "^v^ - p ■ E{x)cos{^) (2) 

Here / = 3 • 10"'^^ gr cm^ is the moment of inertia of rotator equal by orders of magnitude to that 
of a molecule of water; m — 3.2 • 10"^'^ gr is the mass of rotator, close to the mass of a molecule 
of water; p is the dipole moment of a water molecule. We take the electric field of the form 

E = Ax + W (3) 

It is worth nntiner that in our case both t.erms in fhp. kinetic enertrv are « 10"^*^ p.rn. whereas the 



of the potential energy, it is easy to come to the conclusion that Lagrangian function ^ admits 
both finite and infinite regimes of motion. The equations of motion corresponding to Lagrangian 
function ^ read 

/ ^ 

mx = p —— cosUp) (4) 
dx 

lip =— p E{x) sin{ip) 




Figure 2: Motion of a dipole particle represented in the window given by the variables: X 
x-coordinate of the rotator; V velocity of the rotator along x-axis; COS cos{(p) where (p is the 
phase of rotator. Dipole moment and rotator mass p = 1.8D, and 32 • 10^^"' gr. Moment of inertia 
/ = 3 • 10"'^* gram?'. Inserts A, B and C, D correspond to the confined and the infinite motions, 
respectfully. 



Equations (|4]) have only one integral of motion, the energy, and therefore one cannot solve it 
in a finite form, that is by writing down its solution by means of integrals. We have to employ 
numerical simulation for studying it. The key problem is the wise choice of visual representation 
of numerical results. In this respect, it should be noted that equations (j4]) are two equations of 
second order, and therefore their phase space is four dimensional comprising coordinates a;, (p 
and velocities i, ip. We have used variables x, v = x and cos{lp) that provide a kind of three- 
dimensional window on the four-dimensional phase space of equations ([3]). It is important that 
due to the energy integral the values of angular velocity, are bounded, and therefore we can 



the simulation are illustrated in Figl5] We have employed the following values for coordinates and 
fields: 

A initial velocity, phase, and angular velocity: 3-10^ cm/sec, (p — 0, 4-10^"' i/z, respectfully. 
Background field and field gradient W ~ 5 kV/cm and A = 500 kV/cw?. Period of time 
considered 2 • 10""'^^ sec. Amplitude of the particle's oscillation in x several tens of A. 

B initial velocity, phase, and angular velocity: 3-10® cm/sec, = 0, 4-10^'' i/z, respectfully. 
Background field and field gradient W — 5.3 kV/cm and A ~ 500 kV/cm?. Period of time 
considered 10^^° sec. Amplitude of the particle's oscillation in x several tens of A. 

C initial velocity, phase, and angular velocity: 3 • 10^ cm/ sec, = 1.7 rad, Hz, 
Background field and field gradient = 10 kV/cm and ^ = 10 kV/cm? . Infinite motion. 

D initial velocity, phase, and angular velocity: 3 • 10^ cm/sec, Lp — 1.7 rad, 10^^ Hz, 
Background field and field gradient W — 0.1 kV/cm and A = 0.1 kV/cm?. Infinite motion. 

The motion being finite, that is the rotator confined to a finite region of x-axis, depends on 
values of initial coordinates and velocities, and values of the background field, W , and the field 
gradient A. The change of dynamical regimes has the threshold nature, so that small changes of 
fields and initial position may, generally, result in different regimes of motion. 



3 Conclusions 

The numerical analysis of the dynamics of a rotator indicates that there exist regimes corresponding 
to the confinement of rotator to finite regions of phase space. The essential point is that the motion 
takes place in a field that increases linearly as regards the spatial coordinate, when a charged 
particle would move neglecting any boundaries. Thus, the finite motion is due to the interplay 
of translational and rotational degrees of freedom. In this respect it strongly resembles Maxwell's 
pendulum. The latter consists of a fiywhecl and two wires wound round the flywheel axis in the 
same direction and connected with a horizontal support. When released the flywheel goes down 
rotating round its axis due to the attached wires, until it arrives at the lowest point allowed by 
unwinding wires. Then it goes upwards rewinding in the opposite direction. The equations of 
motion are similar to those of the one dimensional rotator, so that one can consider Maxwell's 
pendulum as a kind of mechanical model for the confinement of a particle in channel. 

The model studied in this paper is the qualitative one. It provides the argument in favor of the 
claim that a monomer of water may be confined to a channel, or cavity, inside an ice-like cluster 
of liquid water hydration layers in the vicinity of " membrane water channels" , [7] . Specifically, we 
suggest that the H and H^Q^ ions, which always exist in water, [2], [IHIj could produce a gradient 
electric field that may result in the channeling effect. 

The useful discussion with B.Yu.Sternin is gratefully acknowledged. 
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